The aims of this paper are to show that our single-layer ap- 
proximation for the calculation of the forced oscillations yields 
results that are consistent with the predictions for the synchro- 
nization timescales in circular orbits, T mc ~ a b , thus providing 
a simplified means of computing the energy dissipation rates, 
E. Furthermore, by calibrating our model results to fit the rela- 
tionship between synchronization timescales and orbital sep- 
aration, we are able to constrain the value of the kinemati- 
cal viscosity parameter, v. We compute the values of E for 
a set of 5 M +4 M Q model binary systems with different or- 
bital separations, a, and use these to estimate the synchroniza- 
tion timescales. The resulting T sync f, vs. a relation is compara- 
ble to that of Zahn (1977) for convective envelopes, providing 
a calibration method for the values of v. For the 4+5 M Q bi- 
(^) nary modeled in this paper, v is in the range 0.0015 - 0.0043 
R| /day for orbital periods in the range 2.5 - 25 d. In addi- 
\ tion, E is found to decrease by ~2 orders of magnitude as syn- 
• chronization is approached, implying that binary systems may 
/■""n ' approach synchronization relatively quickly but that it takes a 
much longer timescale to actually attain this condition. The rel- 
| evance of these results is threefold: 1) our model allows an es- 
,— ' • timate for the numerical value of v under arbitrary conditions 
| in the binary system; 2) it can be used to calculate the energy 
■ dissipation rates throughout the orbital cycle for any value of 
C*~) eccentricity and stellar rotational velocity; and 3) it provides 
values of the tangential component of the velocity perturbation 
at any time throughout the orbit and predicts the location on 
the stellar surface where the largest shear instabilities may be 
\Q I occurring. We suggest that one of the possible implication of 
the asymmetric distribution of E over the stellar surface is the 
generation of localized regions of enhanced surface activity. 

Oh 
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Abstract. Stellar oscillations are excited in non-synchronously rotating stars in binary systems due to the tidal forces. 
Tangential components of the tides can drive a shear flow which behaves as a differentially forced rotating structure in a 
stratified outer medium. 



1. Introduction 

A binary system is said to be in equilibrium when the orbit is 
circular, the stellar rotation period and the orbital period are 
synchronized and the spin axes of both stars are aligned and 
perpendicular to the orbital plane. In any other case, oscilla- 
tions are excited and the tidal deformations which are time- 
dependent become inevitably dissipative. We will focus on the 
situation in which the stellar rotation period differs from the 
orbital period, that is, a non-synchronously rotating system. To 
describe the degree of non-synchronicity, we define = cjDqIQ., 
where ojo and Q are the rotational and orbital angular veloci- 
ties, respectively. 

As in the case of the terrestrial tides, the resulting asym- 
metrical shape of the star is due to a second order gravitational 
effect, the differential force. Thus, the rising tidal bulges, oppo- 
site to each other and comparable in size, are associated with 
the so called, equilibrium tide. Additionally, a set of surface 
oscillations commonly known as the dynamical tide, character- 
ized by a wide variety of harmonic frequencies, is established. 
This oscillatory response of the star depends on j3 and on the 
stellar and orbital parameters. In addition, it depends on the ca- 
pability of the stellar material to transfer angular momentum 
and energy between its different layers, which is represented 
by the viscosity, v. 

The tidally induced oscillations have both radial and az- 
imuthal components. If we focus on the azimuthal component 
and keep in mind that the tidal forces decrease as r 3 , with r the 
distance between a mass element and the center of the perturb- 
ing star, it is evident that since an external layer of a perturbed 
star is subjected to a larger tidal force than a layer below it, a 
differentially rotating structure will be produced in any +\ 
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binary system. If v ^0, the interaction between the different 
layers takes the form of a shearing flow, thereby leading to en- 
ergy dissipation. The viscosity and the azimuthal velocity are 
the crucial parameters defining the amount of energy that is 
dissipated, E. 

The long term trend of the angular momentum transfer and 
dissipative effects in binary systems is one in which the rotation 
frequency becomes synchronous with the orbital frequency. 
The theoretical framework for the calculation of tidal synchro- 
nization was developed by Zahn (1966, 1977, 1989), who in- 
voked the tidal friction as a physical mechanism asociated with 
tidally excited oscillations which contribute to the angular mo- 
mentum transfer. In stars with convective envelopes, redistri- 
bution of angular momentum occurs through dissipation pro- 
cesses like turbulent viscosity, retarding the equilibrium tide. 
In stars with radiative envelopes, the mechanism found by Zahn 
(1975, 1977) is the action of radiative damping taking place on 
the dynamical tide. Zahn (1977, 1989) found the characteris- 
tic timescale for synchronization as a function of a, the major 
semi-axis of the orbit; T sync ~ a 6 for stars with convective en- 
velopes and T sync ~ a 8 5 for stars with radiative envelopes. 

On orbital timescales, the tidal oscillations can produce 
observable photometric and absorption line-profile variability 
(see, for example, Smith 1977; Gies & Kullavanijaya 1988; 
Willems & Aerts 2002). Moreno & Koenigsberger (1999) pre- 
sented a simplified model for the calculation of the surface os- 
cillations of a star driven by the tidal interaction in a binary 
system. The calculation is performed for a thin surface layer of 
small surface elements distributed along the equatorial belt of 
the star. The solution of the equations of motion for these sur- 
face elements provides the velocity perturbations in the radial 
and azimuthal directions as a function of time within the orbital 
cycle of the binary. In the non-inertial reference frame used to 
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solve the equations of motion the model includes the gravita- 
tional forces of both stars, gas pressure, centrifugal, coriolis, 
and viscous forces. In the model we compute only the motion 
of the center of mass of each surface element, thus ignoring 
the details of the inner motions within these elements. In this 
approximation, the buoyancy force is not taken into account. 
Finally, we note that the results in this paper are obtained using 
a poly tropic index n=1.5, corresponding to a rigorous convec- 
tive envelope, instead of n= 1 as was used previously (Moreno 
et al. 2005; Koenigsberger et al. 2006), although the difference 
between the results obtained for these two values of n is not 
significant. 



Moreno et al. (2005) extended the model to include the 
computation of the photospheric absorption lines that would 
be produced by the perturbed stellar surface. This ab initio cal- 
culation of the line profile variability due to the forced tidal 
oscillations was applied to the case of e Per, showing that 
from a qualitative standpoint, the predicted line profile vari- 
ability is comparable to the observations. More recently, the 
model was used to analyze the peculiar line-profile variabil- 
ity in the optical counterpart of the X-ray binary 2S01 14+650 
(Koenigsberger et al. 2006), leading to the conclusion that the 
tidal perturbations could be the "trigger" for the periodic ejec- 
tion of wind material at densities that are higher than the av- 
erage wind density. They hypothesize that these effects result 
from the combined action of the shear energy dissipation, the 
radial oscillations and the radiation pressure. The significance 
of this result is that, in addition to providing a tool for studying 
line profile variability in binary systems, the model provides 
a means of constraining the value of the viscosity, which is 
the parameter involved in the energy dissipation mechanisms. 
Thus, although the model is designed to study the variability in 
binary systems on orbital timescales, it is of interest to explore 
whether it yields results that are consistent with the long-term 
behavior, specifically, the tidal evolutionary theories. That is, 
using the value of viscosity and the angular velocity perturba- 
tions obtained from the calculation, an estimate of the rate of 
dissipated mechanical energy, E, can be derived. The timescale 
over which E becomes negligible is a measure of the synchro- 
nization timescale of the binary system. Thus, the results of 
this one-layer calculation of the tidal oscillations may be com- 
pared with the predictions of tidal synchronization theory such 
as those of Zahn (1977, 1989) yielding, on the one hand, a test 
of the Moreno et al. (2005) model and, on the other, a constraint 
on the numerical values of the viscosity. These are two of the 
objectives of this paper. 



The expression for the calculation of the viscous energy 
dissipation rate is derived in Sections 2 and 3; in Section 4 we 
present the results of our model calculations for a 5 M Q + 4 M Q 
binary system at different orbital separations and compare our 
predictions with the synchronization timescales of Zahn; and in 
Section 5 we discuss the implications of this result and present 
the conclusions. 



3 

2. Viscous angular momentum transport in a 
shearing medium 

The approximate treatment of angular momentum transfer be- 
tween the equatorial surface layer and the rigidly-rotating inner 
region of the star is made assuming that the equatorial exter- 
nal layer behaves as a shearing thin disk, in analogy with the 
dissipation mechanisms described for accretion disks (see for 
example, Frank et al. 1985, Lynden-Bell and Pringle, 1974). In 
the assumed differentially rotating medium, tangential stresses 
between adjacent layers are assumed to be the mechanism of 
transport of angular momentum. Included in these is the angu- 
lar momentum transport due to magnetic torques, which has 
been shown to be the dominant mechanism in recent stellar 
model calculations (Petrovic et al. 2005) where the formula- 
tion of Spruit (2002) is used to treat the magnetic field. In our 
calculations, the mechanisms by which transfer of angular mo- 
mentum in the shear flow occurs is englobed in the viscosity. 

We adopt the approach for the transport of angular momen- 
tum that consists of imagining an idealized situation in which 
parts of a ring-shaped outer layer slide over a rigidly-rotating 
inner region. If friction between adjacent layers is assumed to 
exist and if the angular velocity increases or decreases out- 
wards, a net torque will be exerted by the outer ring on the inner 
region. The resulting torques will work towards synchronizing 
the outer and inner region, dissipating energy in the process. 

For this mechanism to be active, the particles in the external 
ring have to be able to interact with those in the inner rotating 
region. However, contrary to the case of many accretion disk 
models where the particle velocity is Keplerian, the velocity of 
different layers on a stellar surface is determined by the stellar 
rotation, the tidal forcing and the viscosity of the material. 

Let us assume that the viscosity in the gaseous shearing 
medium is a turbulent one. Gas elements with small random 
motions at a typical turbulent velocity v tur b, travel a distance 
A-turb (the characteristic length and velocity scales of the turbu- 
lence), in all directions before mixing with the surroundings. 
There exists a typical turbulent viscosity, 

V ~ ^turbVturb (1) 

However, both A tur b and v tm b are unknown parameters. 
A clever solution to this problem was proposed in the early 
1970' s by Shakura and Sunyaev (1973), consisting of a sim- 
ple parametrization of the two turbulence scales. For the length 
scale, they assumed isotropy and characterized it by the typi- 
cal size of the largest turbulent eddies which cannot exceed the 
thickness H of the assumed ring A tur b < H. For the velocity 
scale, a more straightforward scaling could be made. That is, if 
the turbulent velocity were supersonic, emerging shock waves 
would dissipate the energy and reduce the velocity to that of, or 
below, the sound speed. This leads to the choice of a v tur b that 
is subsonic, thus v tur b < C s , and the viscosity can be written; 

v ~ aHC s (2) 

with a < 1 . This is the well-known a-prescription of Shakura 
& Sunyaev (1973). 
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The physical idea of the viscous torque g is to assume that, 
if the gas flow is turbulent to some extent, gas particles of ad- 
jacent layers will be exchanged in the radial direction (Frank 
1985). Since the two different radial flows of material originat- 
ing in both layers will have different specific angular momenta, 
this will cause a net transfer of angular momentum between 
them. 

We can then calculate the magnitude of the viscous torque 
on the equatorial stellar region in terms of a kinematic turbu- 
lent viscosity v, which can be expressed as v = fi/p. Here p 
is the local density of the gas and p. is the coefficient of shear 
viscosity which is related to the component T r<t> of the viscous 
stress tensor in cylindrical coordinates (r, <p, z) by (Landau and 
Lifshitz 1984): 



, 1 du r duj, uj, 

V(-^7 + 1> ) 

r o<p Or r 



(3) 



where u r , u$ are the radial and azimuthal components of the 
velocity field. Now, inserting = aj^r with oj^ as the rotational 
angular velocity of the tidal shear flow, 



~dr~ 



<9(rw^) 



u)</, + r- 



da>s 



(4) 



dr dr 
and neglecting the shear effects in the radial direction, the vis 
cous stress tensor becomes, 

dcOs 



Tr0 = fir- 



dr 



(5) 



This is a tangential force per unit area, and multiplying 
T r by the lever arm r, gives the viscous torque per unit area 
exerted by the sliding ring-shaped outer layer on the rigidly- 
rotating inner region. Integrating over the area of interaction 
(dA = rdcpdz), 

3 du>$ 



g = J rr r<t ,dA = j fir 

In H/2 

= v J r3 ~g^ d( f ) J P dz 



■d(f)dz 



J dr 



d<f> 



(6) 



(7) 



(8) 



where we have used that the density per unit area on the equa- 
torial plane of the equatorial stellar region is 



H/2 

j pd Z 



(9) 



-H/2 



with H the width of the equatorial stellar region. With S = pH, 
the torque can also be expressed as 

2n 

„ da)* 

(10) 

or 

o 

According to this model, no net angular momentum is 
transported unless there is shearing differential rotation; that is, 
the condition for net angular momentum transport is, dco^/dr ^ 
0. 



g = vpHj r'^dcf. 



3. The energy dissipation caused by the 
differential viscous torque 

A natural consequence of the shearing differential rotation is 
the dissipation of energy, E. Let us consider the energy dissi- 
pation caused by the viscous torque acting on the stellar equa- 
torial external layer between r and r + Ar. The net differential 
viscous torque on this layer is 



dg 

g(r + Ar) - g(r) ~ —Ar 



(11) 



Since the torque is acting in the direction of the angular 
velocity aj$, the rate of work exerted onto the external layer is; 



i^-Ar = —{guj^Ar - g—Ar 
or Or or 



(12) 



where the first term of the right-hand side represents the rate of 
transported mechanical energy through the gas by the torque, 
and the term -g(dcj^,/dr)Ar, the corresponding local rate of 
loss of mechanical energy converted into heat in the gas. We 
are interested in the dissipation caused by the viscous torques 
within the gas at a rate g(da>^,/dr)Ar on the external layer of 
width Ar. 

The local dissipation rate on an equatorial surface element 
with an azimuthal width Acp is 



Eeiem = 8~^^ = vpHt 3 I — } ArAcf, 



vp -j r—^— \ rAtpArti 



(13) 



(14) 



Since we require an estimate of the dissipation rate on the 
entire surface of the star, the previous description must be ex- 
tended to include the dissipation produced at polar angles (6). 
An approximate generalization can be performed considering 
that the same dissipation mecanisms apply to a series of rings 
of thickness rA6, situated across the star from the north pole 
to the south pole, and the elements in these rings with an az- 
imuthal angle <p have the same angular velocity as that of the 
equatorial element with the same azimuth <p. This approxima- 
tion was used by Moreno et al. (2005) for the line-profile cal- 
culation. 

From Eq. (14), the local dissipation rate on any surface el- 
ement is 



^ elem 



vp 



<9w . 

r sim 

dr 



r{sinff)A^Ar{rA9) 



(15) 



with a),/, the angular velocity of the equatorial surface element 
in the meridian of the element. Integrating in polar angle, we 
obtain the dissipation rate on a meridional shell, 



- me rid. shell 



2 " 

vpjr^j r 2 A<f>Ar j sin 3 6d9 (16) 



4 f doj,/, , 
3 Vp Y~dr~' ^ 



(17) 
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To obtain the total dissipation rate E, this equation must be 
integrated over the azimuthal coordinate (p. We use the fol- 
lowing approximation: (dco^/dr) ~ (Aw^/Ar),-, where Ato^ ~ 
AV(<f>i) / r i, with r, the radial position of the ;'-th equatorial sur- 
face element, and AV(0,) the difference of azimuthal velocity 
of the element and the inner stellar region. Then, with (Ar), the 
radial thickness of the element, the total energy dissipation can 
be obtained as a sum over the n equatorial surface elements, 




where additionally, r,A0, is simplified by the factor which 
represents the azimuthal length of the element. 

4. Synchronization timescales 

For a super-synchronously rotating star, the transfer of angu- 
lar momentum tends to brake the outer stellar layers, causing 
the total energy of the system to gradually decline. This, in 
turn, drives the system towards synchronization and, therefore, 
|AV(0,-)| -» 0. 

Zahn (1977, 1989) found two different relations between 
the synchronization timescale and the orbital separation in bi- 
nary systems, depending on whether the star's outer envelope is 
radiative, Tzahn ~ («/Ri) 8 ' 5 , or convective, Tzahn ~ (a/Ri) 6 - He 
proposed that for stars with convective envelopes, the most effi- 
cient synchronizing mechanism is due to turbulent viscosity. In 
what follows, we will concentrate on this formulation, a choice 
that is motivated a posteriori by the results we derive and by 
the possibility that, even in stars with radiative envelopes, tur- 
bulence is induced by the differential rotation. 

Zahn (1977) gives in his Eq. 6.1 an approximate character- 
istic synchronization timescale 

Tzahn ~ q~\alR{f ~ 10 4 ((1 + q)/2q) 2 P 4 years (20) 

where q = M2/M1 is the mass ratio of the secondary to the 
primary star and a is the orbital separation in solar units, Ri 
is the equilibrium radius of the star in solar units and P is the 
orbital period in days. 

Adopting our one-layer model for calculating the tidal os- 
cillations, and assuming that all of the energy dissipation asso- 
ciated with the shearing forces occurs on this surface layer, we 
can estimate the timescale required for the rotation period of 
the star to be synchronized with the orbital period of its com- 
panion as follows. If E™' and E^" c are the present and synchro- 
nization stellar rotational energies, computed with the moment 
of inertia of a rigid sphere and the azimuthal velocity at the 
stellar equator, E rot ~(l/2)I(to) 2 ~ (l/5)Mi V 2 ot , then a rough 
estimate of the timescale needed to synchronize the system is 

Ivp^i^iAVicf,,)) 2 } 




j 1 1 1 1 1 1 1 1 1 1 

100 200 

azimuth 



Fig. 1. Top: The difference of azimuthal velocity of each sur- 
face element and the rigidly rotating inner stellar body, AV(<p) 
along the equatorial region for a P or /, =5 d binary system. The 
angle <p is measured with respect to the line connecting the two 
stars. Bottom: The component of the velocity in the radial di- 
rection of each surface element along the equatorial belt for the 
same model. The star is rotating at twice the orbital frequency. 



4.1. Comparison with Zahn's synchronization 
timescales 

In order to compare the above synchronization rates with those 
that are inferred from Zahn's (a/Ri) 6 prediction, our code was 
used to compute values of AV(0,) and r, for a set of binary 
models with different orbital periods. 

The test binary system that was selected consists of a Main 
Sequence intermediate-mass system with Mi =5 M , M2 =4 
M and Ri =3.2 R in a circular orbit. Observational evidence 
(Tassoul 2000; Pan 1997; Mathews & Mathieu 1992; Witte 
& Savonije 1999) provides constraints for the circularization 
timescales of such systems, and T sync is found to be shorter than 
T circ (Claret & Cunha 1997; Zahn 1977; Hut 1981). 

We chose to fix the value of /3 = coo/Q. =2, which means 
that for decreasing orbital separations, the rotational angular 
velocity, too, increases proportionately with the increasing or- 
bital angular velocity £1 Although this might seem like a some- 
what arbitrary choice, it is important to note that it is ft rather 
than V ro( that plays a fundamental role in defining the frequen- 
cies and amplitudes of the oscillations. An intuitive way this 
can be understood is to consider two binary systems with dif- 
ferent a>o, but both with {3 =1. Because (5-1 corresponds to 
the equilibrium condition, no oscillations are present, regard- 
less of how much larger one of the values of too may be with 
respect to the other one. 
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Fig. 2. The azimuthal velocity AV(0) along the equatorial re- 
gion of surface elements with angle <p at a given time, for 
the P or b =2.5 day binary system and with £R=0.01 Ri and 
v =0.0043 Rid" 1 . 

Figure 1 illustrates the example of AV(0) and V r along the 
equatorial belt for the P or j =5 day binary system. The abscissa 
corresponds to the location of the surface element along the 
equator with respect to the line connecting the two stars. The 
largest differences between the surface azimuthal velocity and 
the inner rigid-body rotating region are in general associated 
with the tidal bulges (0 ~0 and 180°) and the locations per- 
pendicular to the axis of the orbit, and it is at these locations 
where the largest energy dissipation occurs. Note that the tan- 
gential component is ~2 orders of magnitude larger than the 
radial component, making it the dominant source of energy dis- 
sipation. 1 When the orbital separation decreases, the azimuthal 
variation is no longer smooth, with shorter scale oscillations 
appearing, as illustrated in Figure 2 where we plot AV(0) for 
the P or b =2.5 day binary system. Because we are considering 
circular orbits, AV(0) shown in Figures 1 and 2 do not change 
significantly as a function of orbital phase and it is possible to 
assign a single AV(0) curve to each of the models. This would 
not be the case if the orbit were eccentric, where strong changes 
occur between periastron and apastron phases. 

The value of the external oscillating layer's density, p, was 
obtained by computing stellar structure models using the code 
described in Yoon & Langer (2005). The value p =10~ 7 g crrT 3 
that we adopted corresponds to the average density in the ex- 
ternal 73 layers of the theoretical 5 M star. These 73 layers 
comprise approximately 6R/R] =0.01, which is the thickness 
of the outer layer that we adopted for the oscillations calcula- 
tion. 

1 It is also the dominant source of photospheric absorption line pro- 
file variability as shown by Moreno et al. 2005. 




10 20 30 40 50 60 70 80 



a (R ) 

Fig. 3. Synchronization timescales (squares) calculated from 
our tidal oscillation model for a fixed f5 =2, and £R=0.01 Ri, 
and for values of v as listed in Table 1 compared with Zahn's 
(1977) relation for stars with convective envelopes (continuous 
line). 



The remaining input parameter that needs to be defined for 
the oscillations code as well as for computing the synchroniza- 
tion times using Eq. (21) is the viscosity, v. Our practice (see 
Moreno et al. 2005) has been to use the smallest value of v al- 
lowed by the code. That is, when the oscillation amplitudes are 
very large, the surface elements may either overlap with or be- 
come detached from their neighboring elements, at which time 
the computation can no longer proceed. This is prevented by 
assigning a sufficiently large value to the viscosity parameter. 
Our initial set of calculations were performed with v =0.005 
R Q day 1 , which is slightly larger than the minimum v allowed 
for the P or b =2.5 day binary system. The results yielded a least- 
means-square fit with T sync Qf(a/Ri) 5 - 7 , which is similar to that 
of Tzahn en (a/Ri) 6 ', but displaced vertically on the logT, vnc vs a 
plane. Thus, the next step consisted in adjusting the values of 
v to make T sync coincide with Tzahn- The resulting values of v 
and energy dissipation rates are listed in the last two columns 
of Table 1, and the synchronization timescales are plotted in 
Figure 3. Note that the 0.0015-0.0043 R Q /d range in v suggests 
that part of the turbulent viscosity is itself produced through the 
tidal perturbations processes, since the closest binaries require 
the largest values of v. 

We thus find that the shear energy dissipation due to the 
tangential velocity components of the surface layer used in our 
approximate treatment of the tidal interaction leads to synchro- 
nization timescales that are consistent with those predicted by 
Zahn's treatment of the star's entire structure for the case in 
which the outer envelope is convective. 
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Table 1. Binary system models for the t swc vs. a calculation. 



p J, 

A orb 


u 


V , 

* orb 


V , 

v rot 


y 


E 


(d) 


(R^) 


(km s" 1 ) 


(km s~') 


(R 2 /d) 




2.5 


16.11 


326.26 


130.00 


0.0043 


1 38x1 34 


5.0 


25.58 


258.96 


64.80 


0.0015 


2.59xl0 32 


7.5 


33.52 


226.22 


43.20 


0.0024 


2.20xl0 31 


10.0 


40.60 


205.53 


32.40 


0.0028 


3.99xl0 30 


12.5 


47.12 


190.80 


25.90 


0.0030 


1.06xl0 30 


15.0 


53.21 


179.55 


21.60 


0.0031 


3.45xl0 29 


17.5 


58.97 


170.55 


18.50 


0.0031 


1.39xl0 29 


20.0 


64.46 


163.13 


16.20 


0.0031 


6.16X10 28 


22.5 


69.72 


156.85 


14.40 


0.0031 


3.02X10 28 


25.0 


74.80 


151.44 


13.00 


0.0031 


1.60xl0 28 



AR=1/100 



1^ = 0.025 



4.2. Dependence on v and 8R/R\ 

The viscosity and the thickness of the oscillating layer are free 
parameters that are not a priori constrained. We find that for a 
given orbital separation, the energy dissipation rates may span a 
range of up to 2 orders of magnitude, depending on the combi- 
nation of v and <5R/Ri . In Figure 4 we illustrate the dependence 
of the energy dissipation rate on viscosity (left panels) and on 
the thickness of the layer (right panels). For P or \, = 2.5 days, 
the energy dissipation rates decrease with increasing viscosity. 
Since the explicit dependence of E on the viscosity is linear, 
(see equation 17) the energy dissipation would be expected to 
increase for increasing viscosity. On the other hand, because 
there is also a strong dependence of the dissipation rates on the 
amplitudes of the oscillations, and larger viscosities tend to re- 
duce the amplitude of oscillation, the net result is a decrease in 
dissipation rates for larger viscosities. We thus see that for the 
closest binary systems, it is this decrease in oscillation ampli- 
tudes with increasing viscosities that dominates the behavior of 
E. 

For the longer orbital periods (P or i, = 10 days), the inter- 
play between the effect of increasing viscosity and decreasing 
oscillation amplitudes becomes more clear. In the middle panel 
of Figure 4 we see how E first increases with increasing v and 
then decreases as the oscillation amplitudes rapidly decline, 
due to the continued increase in v. 

The thickness of the layer also has a strong influence on the 
value of E. In the binary with the shortest orbital period, E in- 
creases systematically by up to 2 orders of magnitude when in- 
creasing 6R from Ri/200 to Ri/50. For longer period binaries, 
however (middle and bottom panels, right), there is a critical 
point at which E reaches a maximum value and then starts to 
decrease with increasing £R/Ri . For the P or \, = 10 days system 
this critical value is 6R C ~ 1/75 Rj while for the P or f, = 25 days 
case the turn around point appears at 5R C ~ 1 / 1 00 Ri. Hence, 
the dependence of dissipation on the thickness of the layer is 
approximately linear up to the thickness <5R C /Ri, with E in- 
creasing monotonically. However, once the layer is thicker than 
<SR f /Ri , the surface layer has more inertia and its azimuthal ve- 
locity amplitudes become smaller and, once again, it is this ef- 
fect that dominates the trend in E. 
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Fig. 4. Energy dissipation rates calculated for different viscosi- 
ties v with a fixed layer thickness (5R=0.01 Rj (left column) and 
for a changing layer thickness with a fixed viscosity v =0.025 
R| d~' (right column). Three binary systems with different or- 
bital periods are illustrated, from top to bottom: P ,b =2.5, 10 
and 25 d. 



5. Discussion 

5.1. The kinematical viscosity 

The energy dissipation rates computed from the results of the 
one-layer stellar oscillation model for binary systems lead to 
the possibility of constraining the value of the viscosity v. This 
can be performed by comparing the value of the synchroniza- 
tion timescale for a given binary system as obtained from our 
model Tgync, with the synchronization timescale curve predicted 
by Zahn's relation for the same system, Tzaim- The difference 
between r !V „ c and Tzaim, if attributed to viscosity alone, allows 
its value to be determined. For the system parameters listed in 
Table 1 we derive v in the range 0.0015-0.0043 R| d" 1 =0.84- 
2.41xl0 14 cm 2 s _1 . If we apply the analogy with accretion disk 
studies where the viscosity is incorporated in the a-parameter, 
an approximate relationship between kinematical viscosity and 
the Shakura & Sunyaev a value is 



kT 1 
v ~ a 



(22) 
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from where 



10 34 w 



FT 



a ~ 0.0981 



(AV„,/km/s)fi(v/Rl/day) 

(r/(ioooo°sr)(/?i/i? ) 



(23) 



Adopting the results illustrated in Figure 2 for the maxi- 
mum azimutal velocity, AV(0), assuming a mean molecular 
weight for a fully ionized gas, ju =0.62, T~10000°K, and using 
v =0.005 R d~', the above equation yields values of a <0.002. 

It is interesting also to compare this result with the values 
of the corresponding sum of diffusion coefficients for rotational 
mixing that are computed by the Binary Evolutionary Code 
(BEC; Petrovic et al. 2005, and references therein). In order to 
make this comparison, BEC was run for a P or /, = 5 d binary sys- 
tem with masses 5+4 M . The computation was stopped after 
an evolutionary time of ~ 10 6 years in order to examine the val- 
ues of these coefficients. The £R=0.01 Ri layer that we adopted 
for the oscillations calculations is equivalent to 73 layers of the 
BEC computation, and the sum of the values of the combined 
diffusion coefficients for rotational mixing 2 from these layers 
is 1.65xl0 14 cm 2 s" 1 (=0.0030 R 2 /d). It is not, however, clear 
how the values obtained for these 73 layers should be com- 
bined in order to compare the result with the v parameter that 
is used in the one-layer calculation, and the sum is most likely 
an overestimate. On the other hand, this value does not include 
the effects due to a magnetic field, which would tend to increase 
the value of the diffusion coefficient for each layer. 



5.2. Evolution of energy dissipation rates 

The energy dissipation rate, which is extracted from the rota- 
tional energy E rot of the star through the action of de-spinning, 
is calculated from the difference between its initial rotational 
energy and the energy at the time of synchronization divided 
by the synchronization timescale. Thus, this represents an av- 
erage of the energy dissipation rate over the time it takes the 
star to reach synchronization. However, the rate of dissipation 
is expected to be a function of time, gradually vanishing as 
the system approaches synchronization. Indeed, Figure 5 illus- 
trates the dependence on /3 of the energy dissipation rate due 
to azimuthal motions, showing that as f3 —>l,E — >0. Figure 
5 corresponds to the model binary system with P or b =5 days, 
v =0.005 R| day -1 , and other parameters as used for Table 1. 

Though we have not addressed the case of excentric bi- 
nary systems in this paper, it is important to note that due to 
the dependence of angular velocity with orbital separation, /? 
is a function of orbital phase. This effect, combined with the 
variation in orbital separation, leads to strongly varying energy 
dissipation rates over the orbital cycle. An exhaustive study 
of the available parameter space is beyond the scope of the 
present paper, but will be the objective of a forthcoming in- 
vestigation. It is, however, worth noting that within the frame- 
work of Zahn's (1977) theory, the synchronization timescales 
are several orders of magnitude smaller than the circulariza- 
tion timescales. Hence, all binary systems should evolve to 
pseudo-synchronization (i.e., average orbital angular velocity= 



10 33 r 

10 32 r 



S ,0- 1 



10 29 t 
10 28 



x x x x x 



1,1,1,1,1,1,1,1,1,1,1 



0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4 

Fig. 5. Energy dissipation rates for different values of /?, for 
a binary system with Mi =5 M Q , M 2 =4 M , Ri =3.2 R , 
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5 d, <5R=0.01 Rj and v =0.005 R 2 d 



2 a-1 



2 see Heger (1998) for a summary and description of the diffusion 
coefficients that are computed 



rotational angular velocity) long before the orbit is circular- 
ized. Thus, large numbers of asynchronous binary systems in 
circular orbits are not expected to exist for stars on the Main 
Sequence. 

5.3. Non-synchronous rotation and surface activity 

An interesting consequence suggested by Figure 5 is that non- 
synchonous binary systems will tend rapidly towards values of 
/3 near unity, but will live for a relatively longer time near the 
equilibrium value, but without having achieved it. Given the 
very large uncertainties in the observable parameters, primarily 
V r0 f sin i and Ri, it is not easy to establish whether a particular 
binary system has actually achieved equilibrium. Thus, there 
may be a large number of non-synchronously rotating binary 
systems in which observable effects due to the tidal oscilla- 
tions may be present. In particular, the peculiar X-ray emission 
that is observed in a small number of B-star binary systems 
in which there is evidence for non-synchronous rotation (Haro 
et al. 2003) may be a manifestation of tidally-induced activity. 
Although the detailed mechanism for converting tidal energy 
into X-ray emission is not known, we have previously specu- 
lated (Moreno et al. 2005; Koenigsberger et al. 2006) that the 
shear produced by the relative azimuthal motions of the outer 
stellar layers may lead to magnetic field generation near the 
stellar surface as well as mass ejection episodes. Thus, the X- 
ray emission could be associated with these processes. The pos- 
sible connection between surface activity and tidal interactions 
has already been investigated in connection with the distribu- 
tion of star spots in RS CVn-type binary systems (Holtzwarth 
& Schiissler (2002). 

Other observational manifestations of this activity may be 
"turbulence" above the expected value for the effective temper- 
ature of the star, and photometric and line profile variability. 
For example, the theoretical absorption line profiles that were 
computed for the optical counterpart of the X-ray binary sys- 
tem 2S01 14+650 required a "macroturbulence" parameter of 
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37.5 km s _1 to match the observations (Koenigsberger et al. 
2006), significantly larger than the thermal speeds expected in 
a B 1 -supergiant star. However, the problem of excessively large 
"macroturbulent" velocities is widely encountered among mas- 
sive stars and, to our knowledge, there has been no systematic 
study made to determine whether the problem is more severe 
in binary stars than in single stars. 



6. Conclusions 

We have shown that a single-layer approximation for the cal- 
culation of tidal oscillations yields results that are consistent 
with the predictions of Zahn's (1966, 1977, 1989) theory for 
the synchronization timescales in circular orbits, T smc ~ a 6 , 
thus providing a simplified means of computing the energy dis- 
sipation rates, E in binary systems. Furthermore, by calibrat- 
ing our model results to fit Zahn's relationship, we are able to 
constrain the value of the kinematical viscosity parameter, v. 
For the binary system models considered in this investigation, 
v -0.00 15-0.0043 R|/day. This is significantly larger than val- 
ues typically obtained from detailed computations of the in- 
terior structure of stellar convective layers, in the absence of 
magnetic fields. 

Our code makes no assumption regarding the energy dissi- 
pation mechanism except that the dissipation processes can be 
described through the parameter v. We show in this paper that 
if we assume that the energy dissipation mechanism is related 
to turbulent viscosity, we arrive at Equation 19, which allows 
us to compute energy dissipation rates using the output of the 
code, and we show that the derived energy dissipation rates are 
consistent with those predicted by theoretical models of stars 
with convective envelopes. However, the code includes no a 
priori assumptions regarding the energy transport mechanism. 
In order to compare the results of the code with the theory for 
stars having outer radiative envelopes, the effects due to radia- 
tive damping need to be incorporated in the derivation of the 
synchronization timescale, which goes beyond the scope of the 
present investigation. 

In summary, however, we have shown that the Moreno et al. 
(2005) one-layer model yields results that are consistent with 
other theoretical calculations of the tidal interactions in binary 
systems. It has the advantage that it can be used to compute 
the tidal oscillations for arbitrary stellar rotational velocities 
and orbital excentricities, allowing an estimate for the numer- 
ical value of v under general conditions in binary systems. In 
addition, the code provides values of the radial and tangential 
components of the oscillation velocities at any time throughout 
the orbit, as well as the location on the stellar surface where the 
largest shear instabilities may be occurring. This allows a com- 
parison with observational quantities, such as absorption line 
profiles and diagnostics of surface activity and should allow a 
methodical exploration of the impact that tidal oscillations may 
have on the surface properties of stars in non-synchronously ro- 
tating binary systems. 
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